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In this paper we study the self-propulsion of a dumbbell micro-robot submerged in a 
viscous fluid. The micro-robot consists of two rigid spherical beads connected by a rod 
or a spring; the rod's/spring's length is changing periodically. The constant density of 
each sphere differs from the density of a fluid, while the whole micro-robot has neutral 
buoyancy. An effective oscillating gravity field is created via rigid-body oscillations of 
the fluid. Our calculations show that the micro-robot undertakes both translational and 
rotational motion. Using an asymptotic procedure containing a two-timing method and 
a distinguished limit, we obtain analytic expressions for the averaged self-propulsion 
velocity and averaged angular velocity The important special case of zero angular velocity 
represents rectilinear self-propulsion with constant velocity. 



1. Introduction 

The study of self-propelling micro-robots is a flourishing modern research topic, striving 
to create a fundamental base for modern applications in medicine and technology, see e. 



Purcell (1977) 


Becker et.al. (2003) 


Najafi & Golestanian (2004) 


Dreyfus et.al. (2005) 




Felderhof (2006)[ Earl et.al. (2007)[ 


Chang et.al. (2007) Alouges et.al. (2008)[|Gilbert at.al. (2010)| 


Golestanian & Ajdari (2008) Alexander et.al. (2009)| Leoni et.al. (2009)| Belovs & Cerbers (2009)| 


Lauga (2011) Romanczuk et.al. (2012) We define self-propulsior 


i as the motion of a 



micro-robot which is subjected to zero external total force. The simplicity of micro- 
robots' geometry represents a major advantage in contrast to the extreme complexity 
of self-swimming microorganisms. This advantage allows us to describe the motions of 
micro-robots in greater depth. The major problem in the designing of a micro-robot is 
the need for an external source of energy to provide for its oscillatory behaviour. Pro- 



posed sources include an oscillating (or rotating) magnetic field (see Dreyfus et.al. (2005) 



Belovs fc Cerbers (2009) Gilbert at.al. (2010) ), an electric field (see|Chang et.al. (2007)), 



and even molecular Brownian forces (see Romanczuk et.al. (2012)[ ). At the same time, 
the major oscillatory forces available in fluid have not been exploited; these are the forces 
caused by fluid oscillations which are imposed by periodically varying boundary condi- 
tions, waves, or turbulence. The ratio of characteristic spatial scales (several microns for 
a micro-robot vs. millimeters, centimeters, or greater scales for flow oscillations) makes 
it clear that the first problem to study is the behaviour of a micro-robot in a fluid that 
oscillates as a rigid body. 

In this paper, we consider the self-propulsion of a two-sphere Buoyancy-dvivcii Dumb- 
bell micro-robot (we call it .RD-robot, see the figure). The whole micro-robot is neu- 
trally buoyant (in order to avoid sedimentation); one of its spherical beads is posi- 
tively buoyant and the other is negatively buoyant. We study two versions of BD- 
robots. In the first one the beads are connected by a rod of prescribed oscillating 




length, in the second one the beads are linked by an elastic spring. First, we study 
the case of a rod and, next, we consider the changes that appear after replacing the 
rod with a spring. A mathematical formulation of the problem leads us to the study 
of creeping motions with time-periodic forces. The problem is solved by employing 
a version of the two-timing method and distinguished limit arguments, developed in 
Vladimirov (2005)[|Vkdimirov (2008))|Vladimirov (2012a)| The approach allows any mo- 



tion of 5£>-robot to be described analytically. Our calculations show that, generally, the 
_&D-robot participates in both translational and rotational motion. Rectilinear transla- 
tional self-propulsion with constant velocity represents a special case of this solution. 
We have calculated the velocity of rectilinear self-propulsion and the ranges of governing 
parameters that correspond to translational motions. 



2. Problem formulation 

The 5f-robot represents a dumbbell configuration, which consists of two homogeneous 
rigid spherical beads of different radii R u , v = 1, 2 connected by a rod of length I (see the 
figure). We study two-dimensional motions of a tree-dimensional dumbbell in cartesian 
coordinates (x,y). The centers of the spheres x^' are described as 

x {1) =X + a, x {2) =X + b, R 1 a + R 2 b = (2.1) 
a = an = riln, b = bn = — r 2 ?n; n = R2HR1 + R2), r 2 = R\/{R\ + R2) 

where X = {X, Y) is the radius- vector of a center of reaction. The axis of symmetry of 
a dumbbell is given by the vector I = x^ 1 ' — x^ , I = \l\. The unit vectors n, cr and the 
angle tp are given by 

f cos ip \ ( -sinp \ . . 

n= . , cr = , er = n v , n^-cr^, n ■ cr = (2.2) 

y simp J y cos<y9 J v v y ' 

where the subscript ip stands for d/dip. The length I is changing periodically 

I = L + sI(t); t = uit: tu = const, e = const (2-3) 

where L is a constant averaged value and I is a 27r-periodic function of r with zero 
average value (throughout the paper a 'tilde' above a function of time denotes that this 
function is oscillating and has zero mean value) . The spheres experience external friction 
forces F M = (f["\f^ v) ) while the rod is so thin (in comparison with either i?„) that 
its interaction with the fluid can be considered negligible. 

We consider the motions of a 5D-robot in a viscous incompressible fluid which, in 
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the absence of the DD-robot, oscillates as a rigid body. These rigid-body oscillations are 
prescribed as a two-dimensional translational spatial displacement £(t) = (£i(t), £2(1")) 
of fluid particles (at infinity in space); the related acceleration is £ tt = uj 2 £ tt , where 
the subscripts stand for related derivatives. The problem can be studied in an oscillating 
(non-incrtial) system of reference, in which a fluid at infinity is in a state of rest. In 
this frame, according to Einstein's principle of equivalence, or according to a related 
transformation of a Lagrangian function, the equations of fluid motion are standard, 
however, they contain an additional oscillating gravity force 

9 = ~^lrr (2.4) 

which causes buoyancy forces f b = —M^g, where the coefficient is equal to the 

difference in the mass of a sphere and the mass of a displaced fluid; M^ v > can be either 
positive or negative. The potential energy of a sphere is IB"' = M^g-x^\ We consider 
a BD-robot of neutral total buoyancy, with total potential energy 

n = n (1) + n (2) = M {1) g ■ i, m = m (1) = -m {2) > o (2.5) 

The problem formulation contains three characteristic lengths: the length of the rod L, 
the radius of the spheres R, and the amplitude of the rod's oscillation a. In addition 
we have the characteristic time-scale T, excess mass M, gravity G, and viscous force F. 
We have chosen these scales as R = (i?i + R%)/2, T = a = sL, F = 6ttt]RL/T, 
g = max|g(r)|, where r\ is the fluid viscosity. The dimensionless variables (marked with 
asterisks) are x = Lx*, t = T t* , Fi = F F* . Three independent small parameters of the 
problem are 

£ = a/L, S = 3R/(AL), m = Mg/F (2.6) 

Below we use only dimensionless variables, however all asterisks are omitted. Note, that 
in the chosen dimensionless units, R\ + R2 = 2. 

We choose the generalised coordinates of _&D-robot to be q — (qi, q 2} 93, 94) = (X, Y, I, <p). 
The motion of the -BD-robot, with a given 1(t) (|2.3[) . is described by the Lagrangian func- 
tion C = C(q, q t ), which includes the constraint (|2.3I) with lagrangian multiplier A 

C{q,q t ) = K-Il + X(l-L-el) (2.7) 

where /C, and II are kinetic energy and potential energy (|2.5p of BD-robot; A represents 
an additional unknown function of time. The Lagrange equations are 

dtdq nt dq n Z^ik^i Qn 

where Q — (Qi, Q2, Q3, Qi) is the generalized external viscous force, exerted by a fluid 
on the .BD-robot. As one can see, we use latin subscripts (i,k = 1,2) for cartesian 
components of vectors and tensors, subscript n = 1,2,3,4 for generalised coordinates, 
and subscripts (or superscripts) fi, v = 1, 2 to identify the spheres. 

We accept, that the fluid flow past BD-robot is described by the Stokes equations, 
where all inertial terms are neglected. In the consistent approximation the mass of a rod 
and masses of the spheres are negligible, hence /C = 0. Therefore (|2.7[) . (|2.8j) . (|2.1|) give 
use to the following system of equations 

F (1) +F (2) =0 (2.9) 
eag er = F a, 2F~ = R 2 F {1) - i? x F (2) (2.10) 
eagn-\ = F~n (2.11) 
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which is supplemented by constraint (|2.3[) . The great advantage of a Lagrangian formal- 
ism is its self-sufficiency. In particular, the conditions of zero force (|2.9[) and the balance 
of torques (|2.10[) appear automatically, while eqn. (|2.11[) allows us to find the reaction of 
constraint A. In (|2.10l) . (|2.11[) we have accepted that m = ea with constant a = 0(1); this 
is our physical assumption, which states that two small parameters e and m (out of three 
parameters in the list (|2.6[) ) are of the same order. Physically, it means that the difference 
between the densities (of each sphere and the fluid) or the amplitude of oscillations of 
the fluid is small (or both these parameters are small). The explicit expressions for 
are 

F« ~ -Rixf ] + 5kS,x1 2 \ ~ -R 2 xf ] + SkSx^ (2.12) 

IS = l 3 Sik = l 2 Sik + Wfci K = R1R2 

Each force F^ v ' represents the first approximation (with the error 0(5 3 )) for the Stokes 
friction force exerted on a sphere moving in a flow field generated by another sphere. To 
construct (|2.12p we use a classical explicit formula for the fluid velocity past a moving 



sphere, see |Lamb (1932)] |Landau fc Lifshitz (1959")] |Moffatt (1996)] Eqns. (f279|) - (f^T2|) 



represent a system of four equations for four unknown functions of time: X , Y, <p, and 
A. For the prescribed I (|2 .3[> . the equation (|2.11j) need not to be considered if we are 
interested only in the motion of the micro-robot and are not calculating of reaction force 
A. For future use, we rewrite (|2.9[) . (|2.10[) as 

X t - 5k§ [X t - R~ i t /4] = (2.13) 
tr-[l t + 5§(R~X t + nl t )] = -2ea(T ■ g/n (2.14) 

where R~ = Ri — R 2 . 



3. Two-timing method and asymptotic procedure 

3.1. Functions and notations 
The following dimensionless notations and definitions are in use: 

(i) s and r denote slow and fast times; subscripts r and s stand for related partial 
derivatives. 

(ii) A dimensionless function, say h = h(s,T), belongs to the class X if h = 0(1) and 
all partial s-, and r-derivatives of h (required for our consideration) are also 0(1). In 
this paper all functions belong to class X, while all small parameters appear as explicit 
multipliers. 

(iii) We consider only functions periodic in t {h £ V : h(s,r) = h(s,T + 27r)}, where 
s-dependence is not specified. Hence, all functions considered below belong to Vf]X. 

(iv) For arbitrary h € V the averaging operation is 

1 pT +2lT 

(h) = — h(s,T)dT = h(s), Vr (3.1) 

2n J to 

(v) The oscillating part of an integral is: 

h T = I h(x 1 s,v) dv f ( f h(x, s, v) dv \ d/i 

Jo 2tt J \J q J 

(vi) The tilde- function (or purely oscillating function) represents a special case of V- 
function with zero average {h ) =0. The bar-function (or mean-function) h = h(s) does 
not depend on r. For any periodic function h a unique decomposition h = h + h is valid. 
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3.2. Asymptotic procedure and successive approximations 

The introduction of a fast time variable r and a slow time variable s represents a crucial 
step in our asymptotic procedure. We choose r = t and s = e 2 t. This choice can be 



justified by the same distinguished limit arguments as in Vladimirov (2012a) Here we 
present this choice without proof, however, its most important part (that this choice 
leads to a valid asymptotic procedure) is exposed and exploited below. We use the chain 
rule 

d/dt = d/dr + e 2 d/ds (3.2) 

and then we accept (temporarily) that r and s represent two independent variables. 
Further more we consider series expansions in the small parameter e and restrict our 
attention to terms which are at most 0(e 2 ). Simultaneously, we keep at most linear in 
S terms. It does not mean that in our setting 5 ~ e 2 , since in all expressions <5 appears 
not separately but as products with various degrees of e. Hence, we do not specify the 
dependence of unknown functions on 5; such dependence reveals itself naturally during 
the calculations. The unknown functions are taken as regular series in e 

X(t, s) = X {t, s, S) + eX^r, s, 5) + e 2 X 2 (r, s, 6) + . . . (3.3) 

with a similar expression for tp. We accept that 

X (s,t,S) = and !po(s,T,S) = while X (s, S) ^ and 6 Q (s, 6) ^ 

which express a basic property of our solutions that long distances of self-swimming and 
large angles of rotation are caused by small oscillations. The application of (|3.2[) to (|3.3j) 
gives 

X t =eX lT + e 2 (x 2T + X Qs ) + 0(e 2 ) (3.4) 

and a similar expression for p. In the calculations below all bar-functions belong to the 
zero approximation, while all tilde-functions belong to the first approximation; therefore 
we omit the related subscripts. 

The successive approximations of equations (|2.13|) . (|2.14p yield: 
Terms O(e ) give the identities = 0. 
Terms O{e 1 S ) lead to 

X T = 0, 1 T ■ W = ~2ag ■ W/k (3.5) 

The use of (|2.2p transforms the second equation to the form p T = — 2ag ■ W/k. Then the 
integration of (|3.5p in the class of periodic functions yields 

X = 0. p = -2ag T -W/n (3.6) 

Terms 0(eS): These terms do not vanish and can be easily calculated. However, as one 
can see below, they do not participate in the leading terms of the average motion that 
appear in the orders 0(e 2 ) and 0(e 2 8). 
Terms 0{e 2 ): Eqns. ([2"T31) i. (|3~ll) give 

X s = -6kR-(&J) (3.7) 

where we have used the integration by parts in the average operation (|3.ip . The use of 
the definition of matrix § (f2TT2|) and jO]) yield (S T T) = 2(Jp T )W. Then (J3J]) takes the 
form 

X a = -SaR" (Tg ■ W)a (3.8) 
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Similarly, from (|2.14j) . we can derive the equation for ip and obtain the system of equations 

X s , = -(iUwfr, 7p s ,=U-yG; (3.9) 
s' = 7s, /j, = 5aR~ , 7 = 2«/k, G = (g[g2) 

U = (T(g ■ a)) = -G x sinTp + G 2 cos Tp, Gi = (7^), G 2 = (Tg 2 ) (3.10) 

where we have used the equality ((g T ■ n)(g ■ W)) = G, which is valid by virtue of 
(|2.1|) . (|3.1[) . and g = (31,52)- One can see, that the dynamics of a dumbbell is determined 
by the values of two parameters /j, 7 and by three correlations G, Gi, and G2. We 
can make a general conclusion, based on p.9j) . that the mean translational velocity is 
always perpendicular to the mean symmetry axes of a dumbbell (X s i is directed along 
the normal vector <r). 



4. Prescribed oscillations of the .RD-robot 

4.1. Unidirectional oscillations of a fluid are not effective 

The simplest motion takes place when G = 0. Physically, it means that fluid oscillations 
are unidirectional, see (|2.4[) . In this case equations (|3.9[) produce the integral 

X-C=-^n or (X-C) 2 = m 2 /7 2 (4.1) 
7 

with a vectorial constant of integration C. This equality shows that X(s') changes along 
a circular path (or along an arc of a circle) of small radius /1/7 = 5kR~ /2. The equation 
for To (|3.9[) can be integrated exactly. For unidirectional oscillations along the y-axis, 
when g = (0,172) (|2.4p . the second equation (|3.9j) takes the form </5 s / = G2COSV3. It can 
be integrated, having an initial value <^(0) = <po, as 

,in-= (l + sm^o)e 2G2S, -(l-sin^ ) 
(1 + sin^ )e 2G ^' + (1 - sin^ ) 

which shows that for s — > 00 we have \lp\ — > tt/2; it means that the axis of symmetry of 
a dumbbell is turning monotonically towards the direction of oscillations. Equation (|4.1[) 
describes the simultaneous change of X along the arc of a small circle. It is clear, that 
in the general case of unidirectional oscillations along any direction (different from y), 
the result is the same: the axis of the dumbbell asymptotically approaches the direction 
parallel to the oscillations, and X changes along the arc of a small circle. Therefore, 
we can conclude that any unidirectional oscillations of the fluid do not result in the 
self-propulsion of the BD-robot. 

4.2. Rectilinear self-propulsion without rotation 

For G 7^ we first consider motions without rotation Tp s , = 0. In these cases the angular 
part of p.9j) gives — G\ sin Tp + G2 cosTp = 7G, which immediately leads to 

Tp = p = — arctan(Gi/G2) + arccos(7G / '"J G\ 2 + G 2 2 ) = const (4.3) 

when | 7 G| < \jGi 2 + G 2 2 (4.4) 

Physically, the restriction (|4.4|) means that the BD-rohot can move without rotation if 
the oscillations I are 'strong enough'. In this case the first eqn. in (|3.9j) gives X s > = 
—j-iGW = const, which shows that the BD-robot moves with constant speed \fJ.G\ in the 
fixed direction cr, which is given by the angle ipo ±ir/2 (|4.3|) . where the sign is determined 
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by the correlation G. Striving for more general results, one can show, that the system 
(|3.9|) can be integrated analytically in the general case ~tp s i ^ 0, with the conclusion, that 
if the parameters satisfy (|4.4[) then a trajectory with any initial data ^(0) asymptotically 
(when s —> oo) approaches the same straight paths (|4.3[) as described above. Exact 
integration outside of the range of parameters (|4.4[) is also accessible analytically; it 
produces motions with rotation \Tp g \ > const, which we do not consider in this paper. 
Let us consider a particular example 

gi = — asinr, <?2 = focosr, I = IqSylit; with a > 0, b > 0, Iq > (4.5) 

7T 

G\ = —al /2, G*2 = 0, G = ab/2, ip = — + arccos(7& / 7o) 

For |76/io| < 1 the 5D-robot propels itself with constant speed \X S >\ = fiab/2 along a 
straight path tp = tpo ± tt/2. It is remarkable, that the self-propulsion speed (iab/2 docs 
not depend on the amplitude ^o- However, one should keep in mind that such solutions 
are available only for 'strong' oscillations, when | Zq I > |7&|; for 'very strong' oscillations, 
when \lo/b\ 00, we have tpo — > it. 



5. Elastic RD-robot 

The above results correspond to the prescribed periodic function 1(t) (|2.3|) . which 
can be chosen arbitrary and represents a given time-dependent constraint. However, in 
practice, the oscillations I produced by the forces exerted from an oscillating fluid on 
the beads, are more interesting. A simple way to consider such oscillations is to replace 
the rod with a spring of stiffness k = const. In this case the dimensionless Lagrangian 
function (|2.7p and potential energy (|2.5p become 

£(<?,<?*) =£-II, n = eag ■ I + k{l - lf/2 (5.1) 

One can check that for this Lagrangian function the equations for total force and torque 
(I2.9p . (|2.10j) remain the same, while the equation for the reaction of constraint (|2.11[) 
must be replaced with 

eag ■ n + kl = —nl t ■ n + 0(eS) (5-2) 

The relation between the problem for an elastic BD-robot and the previous problem 
for a BD-robot with the arbitrary oscillation of a rod is evident: the latter considers all 
possible solutions, while the former corresponds to a special subclass of Z(r) only, that 
appears as the result of spring oscillations. Hence the ability for self-propulsion can only 
'worsen' after the introduction of a spring. The order O(e) equation (|5.2[) produces a 
linear equation for I 

l T + Kl = -7<? • n, K = 2k/ k (5.3) 

It gives us 1(t) which must be substituted into G\ and G2 in (|3.9p instead of an arbitrary 
chosen function I. The rest of the problem remains unchanged. The general solution of 
(|5.3p can be obtained analytically in an integral form, or in the form of a Fourier series. 
Both forms are rather cumbersome and are not considered in this paper. Instead, we 
present an example for a gravity field 7j l = —asinr, 52 = bcosr that coincides with 
(j4~5|) . The related solution of (fO)) is 

7= \-r[(P + KN) sin t + (KP - N) cost] , P = -jbsinlp, N = jacoslp 

1 —I— f\ ^ 
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where an exponentially decreasing complementary solution has been dropped. It leads to 
an explicit formula for U (|3.10j) 



U= 2(1 +K 2 ) ^K(a 2 -b 2 ) S hi(^)-ab ( 5 ' 4 ) 

which determines the system of equations (|3.9[) . In this case we obtain the following 
equations for the motion without rotation (Tp s = 0) 

X s < = -IjfiabW, p = <Po = l; arcsin$ = const, $ = J^^^ + K^j (5.5) 

We can see, that any direction of a rectilinear self-propulsion can be arranged by an 
appropriate choice of a and b. It is also interesting that the speed of self-propulsion, /iab/2, 
does not depend on the spring stiffness k; however, the required condition |$| < 1 shows 
that both a small and high stiffness do not lead to rectilinear motion. Another interesting 
conclusion is: for rectilinear motion to exist, the values of vibrational amplitudes a and 
b can not be chosen close to each other. It means that imposed vibrations (|2.4[) must be 
anisotropic (the circular vibrations with a = b and close to them are excluded). Again, 
for lp s , ^ 0, |$| < 1 the system of equations (|3.9j) . ()5.4j) can be integrated analytically. 
The integration shows that any trajectory asymptotically (for s — > oo) approaches (|5.5[) . 
In the case $ > 1 the full system (|3.9[) . (|5.4[) also allows explicit analytical integration; it 
leads to the motions with rotation \Tp s \ > const, which we do not consider in this paper. 



6. Discussion 



(i) Our choice of slow time s = e t (|3.2p agrees with classical studies of self-propulsion 
for low Reynolds numbers, see Taylor (1951)| Blake (1971) Childress (1981) as well as 



the geometric studies of Shapere & Wilczek (1989) 

(ii) The scale of slow time s = e 2 t implies that in order to obtain physical dimen- 
sionlcss velocities we have to multiply X s i and Tp 8 , (|3.9|) by je 2 . Accordingly, the mean 
translational velocity V = 0(e 2 6), while mean angular velocity SI = 0(e 2 ). 

(iii) We have built an asymptotic procedure with two small parameters: e —> and 
S — > 0. Such a setting usually requires the consideration of different asymptotic paths on 
the plane (e, 5) when, say S = 5(e). In our case we can avoid this additional analysis, 
since 6 does not appear separately, but only in combinations like e 2 S. 

(iv) In this paper we consider only plane motions of a three-dimensional dumbbell. This 
class of motions corresponds to two-dimensional oscillations/gravity (|2.4j) . At the same 
time, for experimental realization, it could be necessary to solve a full three-dimensional 
problem. 

(v) For the first experimental studies of self-propulsion of the BD-robot one can con- 
sider: rigid-body oscillations of a fluid, enclosed within a vibrating container; the viscous 
flows, caused by oscillatory boundary conditions; or the oscillations of a fluid due to an 
external acoustic wave. 

(vi) The calculated velocities of self-propulsion are much smaller than flow oscillations 
and can be even smaller than some secondary flows (like acoustic streaming). In this 
case, one may look for situations when a slow (but permanent) self-propulsion leading to 
new physical effects or providing advantages in applications. 

(vii) It is well known, that an oscillating dumbbell is able to self-swim when an os- 
cillating external torque, exerted on a dumbbell, is present; the related discussion can 
be found in |Felderhof (2006)[ |Felderhof (2007)[ |Friedman (2007)[ The self-swimming of 



a magnetically driven oscillating dumbbell has been studied by Gilbert at.al. (2010) 
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(viii) Our approach is technically different from all previous methods employed in the 
studies of micro-robots. The possibility to describe any motion of a _&D-robot explic- 
itly, shows the strength and analytical simplicity of our method. The studies of different 



micro-robots by the same method (as in this paper) can be found in Vladimirov (2012b) 



Vladimirov (2012c) Vladimirov (2012d) In Vladimirov (2012e) the same method re- 
sulted in a new asymptotic model and a new equation (the acoustic- drift equation) for 
the averaged flows generated by acoustic waves. 



The author is grateful to Profs. M.A. Bees, C.J. Fcwster, P.H. Gaskcll, A.D. Gilbert, 
K.I. Ilin, H.K. Moffatt, T.J. Pcdley, and J.W. Pitchford for useful discussions. 
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